Twisting Electromagnetic Fields with Singular Transformation Optics 
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Transformation optics is a powerful method for manipulating the electromagnetic fields with 
metamaterials which resembles a optical Riemannian space for photons. Analogous to the theory 
of topological defects, we propose a geometric extension of transformation optics by introduction of 
linear topological defect which produces spatial torsion. In order to illustrate the optical effect of 
singular transformation optics, we present an example of application to control the twisting degree of 
freedom of light, the orbital and spin angular momentum of light with an optical cosmic string. Our 
results provide a systematic and versatile framework which may open intriguing venues in angular 
momentum engineering of light with artificial metamaterials. 



Introduction. — In recent years, transformation optics 
(TO) has attracted a great deal of attention because it 
promises a variety of intriguing applications from invisi- 
bility cloaking, electromagnetic wormholes to broadband 
light harvesting SEl • The underlying physical idea 
of TO stems from the formal invariance of the govern- 
ing Maxwell's equations [l[ and Helmholtz equation Q 
under a spatial coordinate transformation. The consti- 
tute parameters of materials e and \i are renormalized 
[sl by the transformation's Jacobian matrix which repre- 
sents the metric tensor of the transformed space. Alter- 
native approaches have also been proposed, based on, e.g. 
space-time transformations [13L direct fields transforma- 
tions [14], and spectra domain|15| transformations which 
generalize the transformation media to include moving, 
bianisotropic and nonreciprocal media. 

So far the rapid growing development has proved TO 
as a powerful and well-established tool for the systematic 
design of artificial metamaterials with specified electro- 
magnetic energy flow distributions, and the same idea 
has made important contributions to advancing the ma- 
nipulation of sound wave [16], matter wave [17], dc elec- 
tric current [l8| and even heat flux [l9|. It is usually 
accepted that the inhomogeneous optical material acts 
as a pseudo-Riemannian manifold, where the propaga- 
tion of photons follows null geodesic determined by local 
spatial curvature @, HqHUI • However, the vector nature 
of electromagnetic waves require some extensions of the 
standard TO to include the twisting degree of freedom of 
light, such as optical spin and orbital angular momentum 
(AM). Horsley [30| has earlier shown that an inhomoge- 
neous isotropic Tellegen medium can be interpreted as a 
Weitzenbock space where the chiral parameter is assigned 
with the contorsion tensor. 

Although the approaches in [3, Ho| seem to capable 
to include chiral metamaterial, no further clues with ge- 
ometric intuitive were given to extend the standard TO 
to embody the twisting degree of freedom of light. More- 
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over, the vast majority of reported TO designs are based 
holonomic coordinate transformation which preserve the 
topology between the virtual and physical spaces. Here 
we explore the extension of TO of Pendry and Leonhardt 
to include spatial torsion and investigate the physical ef- 
fects of torsional media on the propagation of electro- 
magnetic waves, focusing on the generation and physical 
effects of electromagnetic torsion by anholonomic coor- 
dinate transformation. Inspired by the theory of topo- 
logical defects in crystalline solid or gravitation [2j429|, 
we introduced torsion in terms of the spatial gradient 
of the inverse lattice deformation in subwavelength scale 
and naturally arises in Maxwell equations in transformed 
space as magnetoelectric coupling terms. In addition, 
anholonomic coordinate transformation in macroscopic 
scale can transform light beams with planar wave front 
to beams with helical wave front. Numerical fullwave cal- 
culations based on a specific example of linear screw dis- 
locations confirm the theoretical predictions. Our model 
not only extends the geometrical and topological aspects 
of TO but also provides an experimental approach for 
simulating of classical and quantum physics in curved 
torsional space in optical laboratory. 

Tocology and geometry of TO. — Similar to the elastic 
deformation of perfect crystal, method of TO relies on a 
spatial coordinate transformation from a prescribed do- 
main in Cartesian coordinate system to a deformed one. 
The topology of TO will change if the spatial continuity 
is broken under the deformation. Our proposed approach 
is based on the anholonomic coordinate transformation 
in the spatial domain, which is often indicated the emer- 
gence of topological defects. 

As a starting point, we considering a coordinate trans- 
formation from an initial region V in the flat space, 
identified by r = {x a \(x, y, z)} to a transformed region 
V in physical space, identified by prime coordinates 
r' = {x^\{x' ,y' ,z')}. Here the Latin and Roman letter 
are used for the superscript to coincide with the symbol of 
tetrad theory of gravity. The action of coordinate trans- 
formation on point r is generally described by a locally 
special Euclidean group SE{3) [31(, i.e. 
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r' = 0(r)r + T(r) (1) 

where O e 50(3) and T e T(3). Therefore, there are six 
degrees of freedom for coordinate transformation describ- 
ing rotational and translational symmetries in R 3 . From 
the view point of topology, the virtual space and the 
transformed space are generally assumed diffeomorphic 
which means the two space are topologically equivalent. 
The material interpretation of TO regard the domains 

V and V as some elastic objects, which carry systems of 
coordinates together with electric displacement field lines 
in an absolute Euclidean space, the transformed region 

V is just a deformed body with curved coordinated lines 
which is labeled by a metric tensor g^ v (r) = ^abffrrff^ 
and the constitutive parameters of transformed media is 
given by 

= V'diW"'; G^ = e«G a) (2) 

with det(-) denoting the matrix determinant and x, G 
representing (e, ji) and (E, H), respectively. Here and 
below, we use the Einstein summation convention in 
which repeated indices are summed. 




virtual space ^ transformed space *£> ' 



FIG. 1. Schematic illustration of anholomonic coordinate 
transformation in a discrete lattice model. A closed loop 
ABCDA in is mapped to a open curve A1B1C1D1 . The fail- 
ure of the closed loop indicates a line defect along z-axis. 

The key point for going beyond Riemannian geome- 
try comes from the observation that the Jacobian ma- 
trix constructed the tetrad fields in the affine space, i.e. 
e O) = dr/dx^, which related the basis vector in V 
and the basic vector e( a ) in D, 

e M = e M e («) ( 3 ) 

where /i, a = 1, 2, 3 are the spatial indices. Torsion tensor 
is defined in terms of the asymmetry part of the affine 
connection T^ v = (r* — r^ M )/2, where the affine con- 
nection is given by Y^ v = e A • d^e v . Torsion vanishes for 
trivial coordinate transformation which create effective 
Riemannian medium for the photons. By definition, a 
transformed space with torsion demands that the second 



partial derivatives of the coordinate transformation don't 
commutable [29| , that implies the coordinate transforma- 
tion is subject to an anholonimic constraint which breaks 
the translational symmetry, 

d 2 x a d 2 x a 
dx^dx v ~ dx v dx^ ^ ^ ' 

Inequality (j4]) means the violation of Schwartz inte- 
grable condition (29| . From a geometric point of view, the 
constrained coordinate transformation leads to a topol- 
ogy transition form the auxiliary space P to P'. In other 
words, the transformed space V must contained some 
topological defects, such as dislocations in continuum 
limit. 

Unlike the well known meaning of curvature which 
characters how a tangent vector is rotated along a closed 
loop, torsion tensor measures the degree of how the tan- 
gent spaces is twisted about a curve when they are par- 
allel transported, 

dE" = (R^ v E'f + T^)dS^ (5) 

If we define the displacement vector from a perfect lat- 
tice to a deformed lattice by u = r — r'. Encircling the 
defect with an arbitrary closed contour C, the disconti- 
nuity of u is given by 

j (dx ■ V)u = b + x (r - r ) (6) 

where the Frank vector ft represents a rotational action, 
the Burgers vector b denotes the translational action and 
tq the reference point. Figure 1 schematic shows the 
failure of parallelogram around a screw dislocation. The 
emergence of topological defects allows us to distinguish 
different loops connecting arbitrary two points in trans- 
formed space via homotopy classes. For the case of trans- 
formed media with a single linear topological defect, the 
fundamental group is tti(S 1 ) = Z. 

Physical effects of torsion. — Physically, the permittiv- 
ity e and permeability /i, i. e. the metric tensor g, fixes 
the causal structure of the space-time as well as the lo- 
cal optical length; the connections describes the parallel 
transport of free photons. The idea to relate vortices 
with optical dislocations seems intuitively obvious since 
vortices are also known as optical phase dislocations with 
an exjp(U9) azimuthal phase dependence. From TO, the 
electric field in transformed medium is related the 
original E via 

E; = e;E a (7) 

Consider a plane wave in virtual space is transformed 
under the action of e^, the contour integral of E 7 around 
z axis over a closed loop, based on Stokes theorem, can 
be expressed as 

/ dx»% = ( / ds^T^Ea (8) 
Jc Js 
= b a E a 



J5J) indicate that a macroscale topological defect with tor- 
sion can convert a plane wavefront to a dislocated one. 
Specifically, a screw dislocation will leading to a optical 
vortex with topological invariant b. 

To see how torsion is linked with magnetoelectric 
coupling in constitutive relations, we follows the idea 
of standard TO based on formal invariance. Starting 
from Maxwell equations in vacuum, anholonomic coor- 
dinate transformation leading to non-vanishing terms 
^■^v\(d 2 x a / dx v dx\)E a . A simple calculations give out 
the following constitutive relations, 
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with s = i.i= y/det{g)g^ and = ^ vX T p vX . 

Consider the local density of optical chirality C which 
was firstly introduced by Lipkin [32[ |33| 



C = £° E' • V x E' - 
2 

According to (j2j), we have 



B' V x B' 
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C = e^xg^U 



(10) 



(11) 



where U = |(eoE 2 + — H 2 ) is the energy density of light 
field. " /i0 

As an illustrative example, we consider a curved space 
with a topological defect which is described by the fol- 
lowing transformation in cylindrical coordinate [27J 



6' = a0: 
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where a = 1 — 4/i, /3 = b/2ii is the Burgers vector. Obvi- 
ously, this line element turns to be vacuum for a = 1 and 
/3 = 0, a conical space for /? = 0, and a screw dislocation 
or a static chiral cosmic string for a = 1. 

Through this transformation, we can obtain the fol- 
lowing constitutive parameters for the material 



M 

e = fi = — 

cr 



(13) 



with 



M = 



2 + a 2 cos 2 6 sin0cos6[a 2 — 1] 



jsinO 
—jcosO 
a 2 +7 2 



The material parameters (fT3|) describe an optical 
medium with the material distribution breaking trans- 
lational and/or rotational symmetries. To explore the 
optical properties, we perform fullwave numerical simu- 
lations with finite element method (FEM). Figure [2] de- 
picts the electromagnetic field distribution of a Gaussian 
beam incident on a cylinder of defect materials with the 
defect line is aligned with z— axis. It is shown the plane 
wavefront of a gaussian wave is titled in the defect media 
region, and the EM wave intensity in the defect core. As 




FIG. 2. Propagation of a Gaussian beam along an optical line 
defect with (a) f3 = 0.15, a = 1, (b) f3 = -0.1, a = 1 and (c) 
f3 -0.1, a 1.3. 



a result, a gaussian beam is converted to be an optical 
vortex which carrying orbital angular momentum. 

For a specific parameter /3 = 0, coordinate transforma- 
tion (pQ) represents a static cosmic string. In Fig. [3j we 
plot the scattering pattern of a vertically incident gaus- 
sian beam on a cylinder consist of 2D conical media, i. 
e. /i / 1, and /3 = with \i is the deficit angle. For 
ji = 0, the conical medium account for air. When /i de- 
creased, the incident beam is split to two beams with the 
split angle increased with the deficit angle. We noted the 
observed phenomena is an optical analogy of gravitation 
lensing effect by 'cosmic strings'. 




FIG. 3. Scattering of Gaussian beams by cylinders with 2D 
conical geometry with different /x: (a,b) \i = —1.6, (c,d) \i = 
-0.8, (e,f) m = -0.2, (g,h) n = 0, (ij) m = 0.04. 

Conclusion. — We show that the geometry framework 
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of transformation optics, the powerful methodology for 
engineering the electromagnetic fields with metamateri- 
als, can be extended to Riemann-Cartan geometry. Re- 
lying on the anholonomic coordinate transformation, one 
can create optical analogy of topological defect and in- 
troduce torsion tensor in transformation media. In order 
to illustrate the optical effect of singular transformation 
optics, we present an example of application to control- 
ling the angular momentum, both orbital and spin, of 
light with an optical cosmic string. Our results provide 
a systematic and versatile framework which may open 
intriguing venues in angular momentum engineering of 
light with artificial metamaterials. 
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